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Exotic nuclear structures such as bubbles and tori are analyzed through semiclassi- 
cal extended Thomas- Fermi calculations with the Skyrme force SkM*. The variational 
equations for neutron and proton densities are solved fully self-consistently in spherical 
(bubbles) and cylindrical (tori) symmetries. The possible existence of bubble configu- 
rations in some astrophysical scenarios is discussed. The stability of toroidal structures 
against change of quadrupole moment is studied. A global minimum of the energy is 
found in heavy toroidal nuclei. 

1. Introduction 

The idea that stable or metastable niiclei could exist in the form of spherical bub- 
bles was suggested many years agol ^ l ^ l '^ l In the seventies Won^ investigated the 
stability of exotic shapes, namely spherical bubbles and tori, in the framework 
of the liquid drop model (LDM)^ with Strutinsky shell corrections .iS! Wong found 
that within the LDM bubble and toroidal structures were unstable against some 
deformations, but they could be stabilized with the help of shell corrections. He 
addressed his study to known nuclei near the the /3-stability valley. More recently, 
Dietrich and Pomorskl^ have carried out a similar analysis for bubble nuclei with 
large mass numbers. They have found stability islands for spherical bubbles in the 
range of mass numbers A between 450 and 3000. Possible large bubble structures 
have also been analyzed within the mean-field approxi matio n with the Hartree- 
Fock-Bogoliubov (HFB) theory and the DIS Gogny force .'^'^ 

Semiclassical methods of Thomas-Fermi (TF) type are a powerful tool for the 
study of nuclear properties that smoothly vary with the number of particles A. 
Similarly to the LDM the TF methods are free of shell corrections, but they describe 
the nuclear surface in a consistent way.^^ As is known, the self-consistent density 
profiles at the pure TF level end at the classical turning point . This deficiency can 
be cured by using the extended Thomas- Fermi (ETF) methociimSl which adds to 
the pure TF approach fi^ (and fi,^) corrections coming from the non-commutativity 
between the position and momentum operators. The ETF approach allows to obtain 
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self-consistent density profiles with a fall-off similar to the one obtained with the 
mean field Hartree-Fock (HF) or HFB calculations. 

The TF and ETF approaches are free of shell effects and replace the wave 
functions by the proton and neutron local densities as the main variables. This fact 
makes these approaches particularly useful for dealing with calculations with a very 
large number of particles or that involve the continuum (as the description of hot 
nuclei) where complete HF calculations may be extremely complicated. 

Our aim in this paper is to study with some detail the ETF predictions for large 
spherical bubbles and toroidal distributions. It is a step towards a more complete 
description of these structures, that will imply to compute the shell energy with the 
same effective interaction used to calculate the ETF energy. In the second section 
we briefly outline the basic theory. In the third section we present our results for 
spherical bubbles and discuss the existence of two different density distributions 
with almost the same energy per particle for the same nucleus. Some attention is 
paid to the possibility that bubble structures can appear in model astrophysical 
scenarios. In the fourth section we analyze toroidal distributions in deformed ETF 
calculations constrained by a negative mass quadrupole moment. Our conclusions 
and outlook are given in the last section. 



2. Basic theory 

In this work we will use for our semiclassical study of bubble and toroidal structures 
the zero-range SkM * interactiorff^ which is known to give realistic fission barriers in 
actinide nuclei.^m This seems a reasonable choice for describing nuclear structures 
containing several hundreds of nucleons. 

The Skyrme energy density functional can be written as: 

^ — ^kin ^ ^zr ^ ^den ^" ^fin ^ ^so ^" ^Coul: (1) 

where Skin is the kinetic energy including the density-dependent effective mass, 
£zr and Sden are the zero-range central and density-dependent terms which are 
functions of the neutron and proton densities pn and pp only, Sfin depends on Vp„ 
and Vpp and simulates the finite range of the effective nucleon-nucleon interaction, 
and finally £so and Scoui are the contributions to the energy density coming from 
the spin-orbit and Coulomb interactions, respectively. 

In the semiclassical ETF approach up to order, the kinetic and spin-orbit 
energy densities are expressed as functionals of the local density as follows: 

q ^ q Q 
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^- = -|^Ef(W,f, (3) 

where q — n,p for neutrons and protons, respectively. The ratio of the bare nucleon 
mass to the effective mass fq and the spin-orbit potential Wg are the standard 
expressions associated to the Skyrme energy density functional)!^ 
The Coulomb energy density reads 

Seoul = lp{r-}Vcoui{r-} ~ Je2(|)'^'p^/3(f), (4) 

where we obtain the Coulomb potential Vcouiir) by solving the discrete Poisson 
equation AVcouii'r) — --^iTpp{r) using Gaussian elimination. 

From the energy density ([T]) the variational Euler-Lagrange (EL) equations with 
the constraint on given number of neutrons and protons read 

— / dr{£ - ^i^pr, - fipPp) = V— h A— Pq^O. (5) 

Spq J dpq dVpq dApq 

The explicit expressions for the EL equations associated to the ETF energy 
functional up to order can be found in Ref. We solve the set of two coupled 
nonlinear second-order differential equations ^ by means of the so-called imagi- 
nary time-step method (ITSM).[B1 The procedure we follow is described in detail 
in Appendix B of Ref. Our solution of Eqs. ([5]) is fully self-consistent, with 
no initial assumption about the shape of the neutron and proton density profiles 
which are obtained by solving Eqs. ([5]) once two boundary conditions are imposed. 
The EL equations of the ETF model have been solved self-consistently in spheri- 
cal and cylindrical coordinates using the ITSAilSI. The deformed ETF model has 
been succ essfully applied to the semiclassical study of fission barriers with non- 
relat ivist ic ! ^ ^ I ^ ^ I and relativisticff^ interactions. 



3. Spherical bubbles 

In the left panel of Fig. 1 we display the neutron and proton density profiles of a 
nuclear system with N = 560 neutrons and Z = 240 protons as a representative 
example of bubble nuclei. In this case there exist two solutions of the EL equations 
([5|) with very close energy per nucleon. One of the solutions corresponds to a sort of 
hyperheavy nucleus with a re duce d but finite central density, that are called semi- 
bubbles or unsaturated nuclei .l^l^ The other solution, a true bubble, corresponds to 
a distribution where nucleons are concentrated in a relatively thin layer around the 
surface of a sphere. Actually, for a given atomic number Z, there are several values 
of N for which both solutions, semi-bubbles and true bubbles, coexist, similarly 
to the shape coexistence phenomenon in normal nuclei in the actinide region. This 
is illustrated in Table 1 where the energy per nucleon as well as the neutron and 
proton chemical potentials are reported for the Z — 240 isotopic chain from the 
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neutron to the proton drip lines of the bubble configuration. From Table 1 we see 
that semi-bubbles and true bubbles coexist in the range N — 524 (proton drip line 
for semi-bubbles) up to iV=576. Between N — 524 and N — 498 (proton drip line 
for true bubbles) only the true bubble solution is stable against proton emission, 
although solutions corresponding to unstable semi-bubbles (with a positive proton 
chemical potential) can also be found. Above N = 576 and up to iV = 620 (neutron 
drip line for the true bubbles) we only obtain a true bubble solution. 

From Table 1 we also realize that along the Z = 240 isotopic chain the ground- 
state configura tion corresponds to the true bubble solution. It has been found in 
HF calculation^SlHI that the situation is reversed for smaller values of Z, where the 
semi-bubble configuration becomes the ground state. We have checked that this 
also occurs in the ETF calculations. The energies per nucleon in semi-bubbles and 
bubbles reported in Table 1 are much smaller in absolute value than in normal 
nuclei. This is due to the fact that the binding mechanism of these huge nuclear 
clusters stems from a delicate balance between the increase of the surface energy and 
the reduction of the Coulomb energy obtained by putting the protons as separated 
as possible. 

The neutron and proton single-particle potentials (spp) are displayed on the 
right panel of Fig. 1. Roughly, neutron and proton spp follow the corresponding 
densities. The neutron spp vanishes at the center of the nucleus for the true bubble 
solution while it is considerably reduced with respect to its minimal value in the 
semi-bubble case. The proton spp develops strong barriers at the center and at the 
outer part of the true bubbles, whereas the semi-bubble only develops the outer 
Coulomb barrier as expected. 

Heavy-ion collisions at intermediate bombarding energies give strong evidences 
for the formation of highly excited nuclei. After some dissipation due to the pre- 
equilibrium emitted particles and rotation, the nucleus can reach an excited equi- 
librium state with some excitation energy E*. This excited state can be described 
by means of a temperature T if the thermalization time is shorter than the typi- 
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Table 1. The energy per nucleoli in MeV as well as neutron and proton 
chemical potentials in MeV of semi-bubblos and true bubbles for Z = 240 
isotopes. 
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-1.585 


-4.535 


-1.256 


-3.025 


524240 


-4.587 


-1.731 


-0.076 


-4.635 


-1.848 


-1.569 


498240 


-4.673 


-2.428 


1.543 


-4.719 


-2.516 


-0.017 



cal decay scale. This state is metastable with respect to nucleon emission and an 
artificial pressure has to be supplied in numerical calculations for avoiding that nu- 
cleons leave the nucleus. A possible way to do it is to put the nucleus in equilibrium 
with a nucleon gas representing the evaporated nucleons. It is assumed that this 
artificial pressure on the nucleus will not alter the actual physical picture of the 
excited nucleus. As soon as the temperature is not zero, the continuum states start 
to be occupied owing to the thermal Fermi factor. One way of isolating the hot 
nucleus of the continuum states was suggested by Bonche, Levit and Vautherin.IISl 
They found that at a given temperature and chemical potential, there exist two 
solutions of the HF equations. One corresponds to the nucleus in equilibrium with 
its evaporated gas, and the other one to the nucleon gas alone. They proposed to 
study the nucleus by means of a thermodynamical potential calculated as the 
difference of the grand potential associated to each one of these two solutions. 
This subtraction procedure can als o be applied in the semiclassical TF or ETF 
approaches as shown by Suraud.fl^ 

It is worth to point out that when temperature increases shell corrections be- 
come less important and they nearly vanish around T ~ 2 — 3 MeV.^^ This fact 
implies that semiclassical approaches of TF type, which are free of shell effects, are 
very well suited for dealing with calculations at high temperature. On the other 
hand, HF calculations have shown that there exists a limiting temperature 
beyond which normal nuclei become unstable because of the Coulomb interaction. 
When temperature grows, the surface energy decreases faster than the Coulomb 
contribution and, conseq uen tly, nucleons are driven out of the nucleus for a suffi- 
ciently high temperature.^^ 

As far as at zero temperature spherical bubbles become stable against defor- 
mation due to shell effects, it could be expected that spherical bubbles should 
disappear around T ^ 2 — 3 MeV when shell effects are washed out. However, using 
a simple model it has been shown that thermal effects could provide an extra sta- 
bilization against deformation.^^QI We have performed thermal TF calculations with 
the semi- bubble ^''"180 and with the true bubble ^°"274 finding in both cases a lim- 
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Fig. 2. Neutron and proton density profiles at T=0, 2.5 and 5 MeV of bubbles containing 500 
(left) and 900 (right) nucleons. 



iting temperature of about 5 MeV. This limiting temperature in bubbles is smaller 
than the one found in normal nuclei, which lies in the range ~ 8 — 10 MeV. 
This fact can be expected due to the large number of protons in the bubbles con- 
sidered here which favor the Coulomb instability when temperature grows. At TF 
level the excitation energy per nucleon at the limiting temperature is 3.40 MeV/A 
for the semi-bubble ^o^lSO and 2.06 MeV /A for the true bubble 3°°274. This is in 
agreement with the fact that in normal nuclei the limiting excitation energy per 
nucleon decreases with increasing mass number. The fact that huge spherical bub- 
bles cannot sustain a large amount of excitation energy has also been qualitatively 
discussed in Ref. In Fig. 2 the self-consistent TF ne utron and proton density 
profiles corresponding to the liquid plus gas phasJI^G^ at T = 0, 2.5, and 5 MeV 
are displayed. When temperature increases nucleons are pushed from the interior to 
the outer part of the system increasing the rms radius and the surface diffuseness 
of the density distributions. 

The fact that isolated spherical bubbles that are stable against neutron and 
proton emission contain a large number of nucleons makes it difficult to obtain 
such nuclear clusters through heavy-ion collisions. As far as exotic nuclear shapes 

91 

such as rods, pasta, etc. have been conjectured to occur in neutron stars , it 
makes some sense the question whether exotic bubbles as considered here can exist 
in some astrophysical scena rio. To investigate this we make use here of the simple 
approach developed in Ref. ^221 jg assumed that in nuclear matter at subnuclear 
densities the nuclei are located in a lattice which we treat in the spherical Wigner- 
Seitz (WS) approach. We consider a mixture of neutrons, protons, electrons and 
eventually neutrinos in (3 equilibrium inside each isolated WS cell assumed to be 
electrically neutral. In the regime that we are interested in, electrons are extremely 
relativistic and can be assumed to be uniformly distributed in the cell. We solve the 
problem self-consistently at TF level disregarding other particles, shell and pairing 
effects. 

We consider first the outer crust of a neutron star where the WS cell contains 
a nuclear cluster surrounded by a neutron gas with an average density in the range 
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Fig. 3. Neutron and proton density profiles in a Wigner-Seitz cell in the outer crust of a neutron 
star (left) and in conditions prevailing during the gravitational collapse of massive stars (right). 



4.3 X 10^^ g/cm"^ ^ P ^ 2.0 x lO^"* g/cm'^.I^In this scenario temperature is basically 
zero and there are no neutrinos trapped inside the WS cell. Fig. 3 (left panel) 
displays the neutron and proton density profiles obtained with our model assuming 
that the system with Z = 250 is /3-stable and the average density is p = 5.77 x 10"'^ 
fm~'^ in the WS cell. This two conditions allow to determine the radius of the WS 
cell which is i? = 107 fm and the number of nucleons inside the cell: A = 29608. The 
calculated energy per nucleon is 1.846 MeV which is close to the HF ground-state 
energy per nucleon at this average density.^^SI 

A different scenario concerns hot dense matter in conditions prevailing during 
the gravitational collapse of massive stars. One key ingredient for describing the 
collapse is just the EOS. In the relatively low density regime 0.001 fm"'^ ^ P ^ 
0.05 fm~^ the most favored phase consists of nucleons congregated into nuclear 
clusters which can be huge {A < 1000) and immersed in a free nucleon sea. In this 
case the WS cells contain neutrons, protons, electron and trapped neutrinos, again 
under the conditions of charge neutrality and /3-equilibriuni (for details see Ref. 1231 
and references therein). We calculate one possible point of the EOS at a constant 
entropy per nucleon S/A — Ifce, as far as the collapse is fairly isentropic. The WS 
cell contains 855 neutrons and 366 protons corresponding to an average density 
p = 0.05 fm""^ and a proton concentration Y^, = 0.3, which are typical conditions in 
this scenario .1241 The neutron and proton density profiles obtained by minimization 
of the total free energy F (including electrons and neutrinos) are displayed in Fig. 3 
(right panel). We obtain two possible solutions with almost the same free energy 
per particle. One corresponds to a semi-bubble distribution {F/A = 28.86 MeV) 
and the other one to a true bubble embedded in a low density fluid {F/A = 28.73 
MeV). 

From these simple estimates in stellar media we conclude that spherical semi- 
bubbles and true bubbles could in principle exist in the considered astrophysical 
context. Whether these bubbles are stable against deformations demands further 
investigations. 
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4. Toroidal structures 

Long-lived nuclear structures beyond the island of stability may be reached if other 
topologies different from the sphere are considered. Simulations using Boltzmann- 
Uehling-Uhlenbeck (BUU) transport equations have shown that exotic shapes such 
as bubbles and tori could be created in central heavy-ion collisions!^ Some ef- 
fort has been addressed to find experimental evidences of the breakup of toroidal 
distributions produced in heavy-ion collisions 

Nuclei of toroidal shape were first analyzed by Wong in the framework of the 
LDM plus Strutinsky shell corrections!^ It was found that there exists a minimum 
of energy for such a type of shape which, however, is usually unstable against multi- 
fragmentation. Within mean-field theory using Skyrme fo rce s the toroidal structure 
in superheavy nuclei was predicted by Nazarewicz el al.^^H More recently, axially 
symmetric toroidal superheavy nuclei with Z > 130 have been analyzed by one of 
us within the HFB theory with the Gogny DIS force.^^ A global minimum of the 
binding energy in toroidal distributions is obtained by performing constrained HFB 
calculations with a relatively large and negative mass quadrupole moment Q2 ■ It is 
also found that for mass numbers A large enough, the superheavy toroidal nucleus 
becomes the ground state and it can be as much as 200 MeV more bound than the 
corresponding spherical nucleus. 

We perform here deformed ETF calculations including ^^-corrections with 
the SkM* force to analyze semiclassically toroidal nuclei. The Weizsacker term 
^ /g(Vpg)^//9g of thc kiuctic energy density given in Eq. ([2]) is multiplied by a 
factor 1.6 in such a way that the ETF surface energy calculated with the SkM* 
interaction equals the HF value. " 

Fig. 4 displays the potential energy surface (PES) as a function of the quadrupole 
mass moment Q2 for systems with = "^^^122, ^20132, 352^^44 ^-^^ 4i6^g4^ g^^^j^ 
systems would hypothetically be created in ^o^Pb -I- ^°Zr, ^ospj^ _|_ nssn^ ^ospi^, _|_ 
^'^''Sm and ^^^Vh -\- ^''^Pb collisions of stable nuclei, respectively. In all the analyzed 
cases we see that the PES consists of two branches that cross at a given value of 
Q2- The crossing point is shifted to more negative values of Q2 when the mass 
number of the system increases. The dashed branches in the Fig 4. correspond to 
a compact oblate distribution and the solid ones to a toroidal shape. These two 
types of nuclear distributions can be seen in Fig. 5 where 2-dimensional plots of the 
self-consistent density of the nucleus ^^^164 are displayed for the compact solution 
at Q2 — —250 b and for the toroidal structure at quadrupole moments Q2 = —250 
b, —650 b and —990 b (ground state). 

For the lightest system 29^122, the ground state corresponds to a compact oblate 
shape with comparatively small negative deformation. However, when A and Z 
increase, the minimum point of the toroidal branch decreases and it becomes the 
ground state for the ^20132 system, as well as for the heavier systems. From Fig. 4 
one also observes that when the mass and the charge of the system increase the 
ground-state energy is more negative as compared with the energy of the spherical 
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solution {Q2 = 0). For heavy systems the minimum becomes more flat and its 
position is shifted to more negative values of the quadrupole moment making the 
toroidal nucleus unstable against multifragmentation. The PESs of toroidal nuclei 
calculated in the Skyrme ETF method (without shell corrections) are in qualitative 
agreement with previous findings in the mean-field HFB theory with the Gogny 
force, althou gh t here are differences in the absolute values of the binding energies 
of the system!^ 

In a small range of quadrupole momenta around the crossing Q2 values we find 
numerically two possible solutions of the EL equations ([5|) for the same value of 
Q2. These two solutions have close binding energies. One of them belongs to the 
compact branch and the other one to the toroidal branch in Fig. 4. (This situation 
resembles the one discussed previously for spherical bubble distributions, where it 
is possible to find simultaneously semi-bubbles and true bubbles for some values of 
mass and atomic numbers.) In the case of the ^^^164 system, the two solutions are 
foimd in the range —350 ^ Q2 —200 b (see Fig. 4). Two-dimensional contour 
plots of the compact and toroidal density distributions at Q2 = —250 b for ^^^164 
are depicted in the upper panels of Fig. 5. 

It is interesting to inspect the modification of the calculated PES induced by 
a finite temperature T and by the rotation of the torus around a radial axis with 
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Fig. 5. Two-dimensional density plots in the xz plane for the nucleus '*^®164 at several mass 
quadrupole moments Q2- The configurations are axially symmetric around the z direction. 



angular momentum Lj. . The stability of hot rotating toroidal structures was studied 
in the past by Wong^in the framework of the LDM. In the left panel of Fig. 6 we 
display our ETF resuhs for the case of the ^^^16A system at T = and T = 1 MeV, 
and when the system is rotating with angular momentum = 150 h. To simplify 
the figure, we use the same line style for the compact and toroidal branches of the 
PES in each studied case and do not draw the points in the crossing of the two 
branches. 

One observes that for small values of T the excitation energy of the '^^^164 
nucleus is roughly independent of the deformation Q2. With a temperature T = 
1 MeV the energy is shifted upwards about 55 — 60 MeV along the whole PES. This 
value can be seen as the difference between the dashed and solid curves in the left 
pan el of Fig. 6. We treat rotation in the rigid-body approximation as described in 
Ref. The excitation energy gained by the ^^^164 system under rotation is more 
important for more compact (i.e. less deformed) distributions, as it can be expected 
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Fig. 6. Thermal and rotational (left) and surface energy (right) effects on the PES of the ■*1^164 
nucleus as a function of the mass quadrupolc moment Q2- 



from the reduction of the rigid-body moment of inertia when Q2 increases (i.e. is 
less negative). In the present example, the energy gain of the system with angular 
momentum — 150 fi. with respect to the non-rotating system is around 30 — 35 
MeV in the compact branch of the PES, while it decreases from about 20 MeV at 
Q2 — —350 b to only 7 MeV in the region around the minimum of the PES. 

It is to be noted that the PES of toroidal nuclei depends on the value of the 
surface energy associated with the nuclear interaction. To exemplify this, we display 
in the right panel of Fig. 6 the PES of the ''^^164 system obtained multiplying the 
Weizsacker term in Eq. ([2]) by a factor of 2 (to reproduce the LDM fission barrier 
of ^■*°Pu|I^ instead of by the 1.6 factor used in the previous calculations. It can 
be seen that in this case the minimum exhibited by the PES is more marked and 
shifted towards more compact (less negative Q2) configurations. Consequently, it 
is expected that nuclear forces with large surface energy make toroidal nuclei more 
stable against multifragmcntation. 

5. Summary and outlook 

In this work exotic nuclear shapes, namely bubbles and tori, containing several 
hundreds of nucleons are analyzed by means of the semiclassical extended Thomas- 
Fermi method and the Skyrme force SkM*. Calculations are performed by solving 
fully self-consistently the variational Euler-Lagrange equations in spherical and ax- 
ial symmetries. The ETF method incorporates in a consistent way the description 
of the nuclear surface and is free of shell effects. Although shell effects are crucial 
for the stability against deformation of these huge and exotic nuclear clusters, im- 
portant information about the smooth variation with the number of particles of 
some global properties can be obtained from the semiclassical analysis. 

As it happens in the HE calculations, for some neutron and proton numbers two 
different solutions of the ETF equations, namely semi-bubbles and true bubbles, 
are obtained. For a large enough number of nucleons the true bubbles are the 
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structure with minimal energy. We have also explored thermal effects in spherical 
bubbles. In this picture, spherical bubbles can exist at temperatures above the one 
for which shell effects disappear but with a limiting temperature smaller than for 
normal nuclei. Using a simplified model we have also explored the possible existence 
of spherical bubbles in astrophysical scenarios. From a semiclassical point of view 
it is found that semi-bubbles and true bubbles could exist in the outer crust of 
neutron stars and in hot and dense stellar matter in conditions prevailing during 
the gravitational collapse of massive stars. 

We have performed deformed Extended Thomas- Fermi calculations for toroidal 
distributions of nucleons. Threshold values of the mass and atomic numbers exist 
beyond which toroidal nuclei become the ground state of the system. With increas- 
ing mass and charge the ground state of toroidal nuclei is more deformed and more 
bound with respect to the spherical solution, undergoing ultimately decay by mul- 
tifragmentaion. Effects of temperature, rotation and value of the surface energy on 
the PES of tori have also been discussed. 

The semiclassical study presented here is a step towards a consistent description 
of bubbles and tori which has to include shell corrections computed from the semi- 
classical single-particle potentials. In this way it is expected to obtain a quantal 
description of these exotic nuclear clusters that will allow to know separately the 
smoothly varying and the oscillatory parts of some nuclear properties. The prelimi- 
nary results obtained in astrophysical scenarios are very encouraging. An important 
point is to investigate the stability against deformation of bubbles embedded in a 
gas containing nucleons, electrons and neutrinos. 
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